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Abstract 

We propose an analytic method for predicting the large angle CMBR 
temperature fluctuations induced by model textures. The model makes 
use of only a small number of phenomenological parameters which ought 
to be measured from simple simulations. We derive semi-analytically the 
C'-spectrum for 2 < / < 30 together with its associated non-Gaussian 
cosmic variance error bars. A slightly tilted spectrum with an extra sup- 
pression at low / is found, and we investigate the dependence of the tilt 
on the parameters of the model. We also produce a prediction for the two 
point correlation function. We find a high level of cosmic confusion be- 
tween texture scenarios and standard inflationary theories in any of these 
quantities. However, we discover that a distinctive non-Gaussian signal 
ought to be expected at low /, reflecting the prominent effect of the last 
texture in these multipoles. 
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1 Introduction 



In the recent past, a large amount of work has been directed towards predicting 
the CMBR temperature anisotropics associated with inflationary scenarios (see 
[0 for a review), and with the various topological defect scenarios (e.g. 0, 0). 
Inflationary scenarios are far better explored in this respect. The main reason 
for this is that it is easier in these models to plead ignorance of the detailed 
mechanism responsible for the density fluctuations of the Universe. Quantum 
fluctuations in the metric are a fuzzy subject, and hand-waving arguments for 
Gaussian fluctuations with a particular type of spectrum are normally considered 
acceptable. Once this is done, one is free to take advantage of the well developed 
and sophisticated industry of methods aimed at treating Gaussian fluctuations 
[0. A precise prediction for the outcome of concrete experiments follows easily, 
making inflation popular with experimentalists. The state of affairs in topological 
defect scenarios is rather different. No lack of fundamental physics understanding 
hampers working out in rigour the defect network evolution and the fluctuations 
they induce on the matter and radiation of the Universe. The price to pay for 
our honesty is that a serious treatment of the problem is computationally over- 
powering. A reliable link between the defect network and the perturbations it 
induces is undoubtedly still missing. To top off the trouble, it turns out that 
the fluctuations induced by defects are often non-Gaussian. Setting up a data- 
analysis framework geared towards non-Gaussian fluctuations is virgin ground 
(but see |§, ^ |, |^ for Gaussianity tests). Overall, even though defect scenar- 
ios have so far been poor in solid predictions, they have a potentially polemic 
value not to be discarded. Non-Gaussian data analysis does not exist partly be- 
cause non-Gaussian theories are too vast a class to be treated in one go. Defect 
scenarios provide a motivated and topical example with which one can rehearse 
non-Gaussian data analysis techniques and experimental strategies. 

The main technical complication in predicting CMBR fluctuations in defect 
scenarios arises from the fact that ^ maps photograph the whole sky (that is, our 
past light cone clipped by the last scattering surface) and not just its intersection 
with the last scattering surface. Hence the defect network, surrounding matter, 
and CMBR photons have to be evolved through a large number of expansion 
times. The problem is usually tackled by means of full-blown simulations 0, |^. 
These require access to supercomputers, and even so face severe limitations. Im- 
portant issues like cosmic variance in defect scenarios have hardly been addressed. 
A possible computational short cut can be found if one is to assume scaling. The 
defects network is expected to scale, by which one means that it looks statis- 
tically the same at any time, once the network length scale is equated to the 
horizon distance. This suggests dividing the sky cone into cells corresponding to 
expansion times and horizon volumes. Scaling makes it plausible that the CMBR 
pattern in each of these cells is statistically the same, once some angular scaling 
procedure is applied to it. This being the case, all that is required from the sim- 
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ulations is good statistics for the CMBR spots induced by defects inside horizon 
size boxes during one expansion time. Such simulations are considerably easier 
to perform. A prototype of this type of study is 0. In the case of textures, the 
short-cut which we have described has been adopted by where a "scaling-spot- 
throwing" process was implemented on a computer. The spots were all derived 
from model SSSS (self-similar and spherically symmetric) collapses. In this pa- 
per we advocate the use of analytical techniques for scaling-spot-throwing. An 
analytical approach allows the derivation of exact formulae for the C', C{6), and 
their cosmic variance, from which interesting general properties can be derived. 
These properties are naturally less evident in computer simulations, unless one 
explicitly sets out to find them. Furthermore, an analytical approach allows the 
treatment of a more general case than 0. The method proposed here applies to 
quasi-circular spots with any profile. The scaling size (ps), number density (n), 
and brightness factor (a) of the spots are also left as free parameters. We should 
admit that only the brightness factor is a truly free parameter in the texture 
model. All others should be considered as phenomenological in origin, eventually 
to be fixed by simulations. As they still constitute a controversial matter, we 
choose to consider a broader class of model spots. We can then sort out which 
spot properties it is really important to measure from the simulations. 

The plan of this paper is as follows. In Sec.|^ we start by deriving the statistical 
properties of scaling texture model collapses as they appear in the sky cone, and 
from these a formal solution for the joint distribution function of the CMBR 
is given. Using this solution, in Secs. ^.l] and we derive exact formulae for 
the and cr(C;) in texture models. We prove the general result that the cosmic 
variance in the Ci is larger in texture scenarios than in Gaussian theories. The 
excess variance is negligible for high /, but significant at low /. In Secs.|3.2| and 
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a large class of spot types. We find that the spots' profile, intensity, and number 
density affect only the spectrum normalization. The spectrum shape is largely 
controlled by the spots' scaling size Ps alone. In general, the spectra are slightly 
tilted (and we numerically compute the dependence of the tilt on p^), but with 
an extra suppression of power at low /. In Sec.^ we derive expressions for the 
two-point correlation function C{6) and its cosmic variance. Although closer to 
experiment, C{6) appears to be a very bad discriminator between texture and 
inflationary scenarios. The most original result proved in this paper is presented 
in Sec.^, where we ask how many textures are responsible for a given a|^. For 
high /, this number goes like I, suggesting Gaussian behaviour. For low / we find 
that the are mostly due to the effects of a single texture: the last (that is, 
the closest) texture. This ought to impart a peculiar non-Gaussian signal to the 
low a^, a fact which we prove explicitly. The purpose of the sequel to this paper 
[ pUd is to devise a low-/ data- analysis technique (in the moulds of |TI|]) capable of 
fully characterising this effect. 
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2 Analytical spot-throwing 



In topological defect scenarios the CMBR fluctuations at large angles are due 
only to the ISW (integrated Sachs- Wolfe) effect caused by the defect network 
(0). For a standard recombination scenario the validity range for this approxi- 
mation is 2 < / < 30. If early reionization occurs, the upper limit for this range 
becomes the value of / corresponding to the angular scale of the horizon at "last 
scattering" (defined in some conservative way). As it is not clear which is the 
most sensible scenario we take as a working hypothesis the choice which maxi- 
mizes the applicability of our model: standard recombination. Now if we assume 
that the statistical scaling of the texture fields extends to their associated time- 
dependent metric perturbations, then we may expect the ISW spots themselves 
to scale. This allows the direct modelling of the statistics for texture induced 
spots as they appear in the sky. From these, we can then write an expression for 
the joint distribution function for 2 < / < 30. For simplicity, we will assume 
that the texture spots are independent (not required by the formalism, but a 
good approximation). Then each texture angular position in the sky is uniformly 
distributed. Furthermore its time position in the sky cone can be associated with 
a modified Poisson process in time, in which the probability of an event is allowed 
to vary. Let n be the average number of texture spot producing configurations 
per horizon volume. Then their average volume density is 



and using the relation between the proper area dS^ the time of emission t, and 
the apparent solid angle dVt: 



we find that the angular density of spots in the sky per unit of emission time is 



It is convenient to introduce the variable y = \og2(to/t), the expansion time 
number as we go backwards in time. Each Ay = 1 represents one expansion 
time, starting from yo = (here and now) and going as far as yis = log2(to/^/s)5 
when we hit the last scattering surface. In the standard recombination scenario 
yis ^ 15. In terms of y the spot density in the sky is 






(3) 




(4) 
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Let us now introduce a Poisson-type process in y, with a variable probability 
density N{y). If Pnijj) is the probability of having exactly n textures at y' < y, 
then 

^ = iV(l/)(Pn-l-p„) , P-1=0. (5) 

This hierarchy of equations fully specifies the statistics of the textures' time- 
position in the sky. From the Pn{y) one can also derive the probability density of 
the n^^ texture position ?/„ by means of 

PniVn) = Pn-l{yn)N{yn) ■ (6) 

However one should be wary of the ?/„ variables, as they are not independent. Not 
only are their ranges constrained hj yi < y2 < ■■■ < yn but also their statistical 
dependence is patent in the conditional distribution functions. From and (|^) 
one has that 

Pi(yi) = iV(yi)e-r^^^(^) (7) 

but then 

P2{yM = N{y,)e"^:'''''''^ (8) 

which is dependent on yi. From these expressions we can derive the joint distri- 
bution function 

^12(1/1,2/2) = N{y,)N{y,)e-Io"''y''^^^ (9) 
which misleadingly factorizes. More generally one has 

Pi...n{yuy2,...,yn) = N{y,)N{y,)...N{y^)e-^T^y^^y) . (10) 

The statistics contained in P„ or p„ refer to a sequence of textures ordered in y. 
Ordering will be crucial in Section ^ when we uncover the last texture. However 
for some other purposes (like the derivations in Section ^) we may simply consider 
unordered textures. Then we have a total of Nt = Jq''" dy N{y) ^ 1 textures in 
the sky. Their positions ?/„ are independent random variables with a distribution 

PM = ^ (11) 

for all n, which is considerably simpler. 

We now model the texture spot patterns. These ought to follow some sort of 
statistical scaling law, where the scale angle is 

r(y) = arcsin (min (l, ^^^)) (12) 

for a pattern laid at time y. Here ps is the scaled impact parameter {ps = p/ (3t)) 
giving (Ts = npl, the scaled cross section for photon anisotropics. Causality 
requires that < 1/2. It may happen that ps is a random variable itself, but 
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we shall ignore this complication. For small, nearly circular spots with profile 
W^{6,y), scaling implies that 

W'{e,y) = W{x) (13) 

with X = 9/9^ . This is a sensible scaling law even for y < 1, although for y < 1 
the exact scaling law should be more complicated, as the photons propagating 
through the scaling metric no longer move along parallel trajectories. In partic- 
ular, a y < yo = 31og(l +Ps), we hve inside a texture and are enclosed by its 
pattern. Unless {2nps)^ ^ 1 textures like these are very rare, and the approxi- 
mation (1^) should produce good enough final results. 

We seek to write the aj^ coefficients for skies filled with nearly circular spots 
described by the above statistics. For each of these spots let us first point the 
2;-axis at their centre, and perform the spherical harmonic decomposition in that 
frame. We obtain an expression of the form 



d'M = anW'\yn).^'^^^.S^mo + 4 (14) 

where W^^ is the Legendre transform of the scaled profile: 

W'^ = 27r dz P\z)W'{z, Vn) (15) 

and is a perturbation induced by the eccentricity. This perturbation is negli- 
gible, even if the spots are not very nearly circular, but / ^ [1/6** (?/„)]. The extra 
factor a„ in (plSj) accounts for the texture brightness, and is a random variable 
with a 7/n-independent distribution function. Performing a rotation to a general 
frame, where the n^^ texture has coordinates and summing over all textures 
one finally obtains: 

al = Y.^nW^\yn)y^{^n) . (16) 

n 

The variables a„ are independent but equally distributed. Their distribution has 
to be determined from the simulations. The VLn are uniformly distributed, and 
from them one can determine the y^(fi„) distributions (but with care, as the 
are not independent). Finally from (p^OD, (|12[), and ( p^D one can determine the 
distributions. The W^^ makes the various terms in (|16D dependent variables, 
rendering an analytic approach for the distributions through (|TBp unfeasible. 
In any case finding the marginal distributions of the would not be the end 
of the story, since the o}^ are necessarily dependent random variables (see |]TT|). 
Hence we still would have to find the joint distribution F(a^j, a^^^ •••) ^ com- 
plete solution to the problem of cosmic variance in texture scenarios. Expression 
(|T6|) constitutes a formal solution to this problem which, despite all the obstacles 
to an analytical solution, is particularly suited to Monte Carlo simulations. It 
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separates the individual effect of each texture on each of the al^. It also factor- 
izes each texture contribution into three factors. The are purely geometrical 
factors which should be present in any SSS (statistically spherically symmetric) 
theory (see for a definition). The W^'' are a reflection of the structure of the 
sky (foliation into expansion times and horizon volumes) and should be present 
whenever there is scaling. Finally, the a^s factor out all that is peculiar to the 
texture model under consideration, and will have to be measured from detailed 
texture simulations. 



3 The angular power spectrum 
3.1 An expression for the 

Throughout this paper we will use the notation Ci for the angular power spectrum 
of realizations, and C' for its ensemble average. Whereas Ci is a random variable, 
is a number. In any SSS theory one has 

< >obs= C^^\i5^, . (17) 

The can be estimated from the observable 



C, 



1 ' 



I E (18) 



2/ + . 

m=—l 



since < Ci >obs= in any SSS theory. Cosmic variance introduces error bars of 
size (j{Ci) into this estimate. We start by deriving an exact expression for the 
From (|T6|) and (ITSl) we have that 



=< Ci >= — ^ 
21 + 

but since 

we have simply 



tEE < X wiwi, X Fjr(^]„)ll(^]„o > (i9) 

m nn' 



< y^*i^n)Yi{a.,) >obs= ^ (20) 



& = ^I^E <WL'> . (21) 



Computing the averages using the unordered statistical description (|TTD we can 



write 

,2 



C' = — / dyN{y)W'^\y). (22) 



47r 7o 

If Ps is a random variable with distribution f{ps), then expression (^2|) still holds 
but now with W^"^ replaced by / dpsf{ps)W^'''^{y,ps). Expression (|2^) lends itself 
to be written as 

&= r'dyC\y) (23) 
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Figure 1: The density C^(y) multiphed by /(/ + 1) for / from 2 to 30 for (n = 
l,Ps = 0.4) and (n = l,Ps = 0.08) textures with Gaussian profile spots. 



thus defining a density 



measuring the contribution of textures living at time y to a given CK In Fig. |1] 
we have plotted /(/ + 1)C'(|/) for (n = 1,^^ = 0.4) and (n = l,ps = 0.08) textures 
with Gaussian profile spots, for / from 2 to 30. These plots are to be confronted 
with the common belief that a given is mostly due to textures living at a time 
when the apparent size of the horizon equals the multipole angular scale. We 
find instead that C\y) peaks for y: Q^{y) ~ J(j+T)^ that is, when the textures' 
apparent size fits the multipole angular scale. Furthermore the peaks comprise 
several generations and an accurate prediction for the C' spectrum can never 



neglect any of the integrand in (23). Changing n only affects the normalization 
of C\y). Reducing ps not only decreases the normalization but also shifts the 
peaks to the left necessarily affecting the shape of the integrated CK In general 
as / increases the peak's heights (multiplied by /(/ + 1)) are slightly reduced, but 
their width increases. The delicate balance of these two effects will determine 
whether the C' spectrum is flat or not. 



3.2 Profiles and reduced a and ps 

It is not clear what proflle W real texture spots have. Simulations [0, |] suggest 
that the well-known SSSS proflle may be oversimplifled. We flnd that the exact 
proflle form has little impact on the C'-spectra. For intermediate /, proflles do 
not matter at all, and for low I they induce differences easily confused by cosmic 
variance. The issue of matching free parameters for different proflles is, however, 
non-trivial. We have considered the following proflles: 
l.-Hat profile. Deflned by 

W\x)=H{l-x) (25) 

where H is the Heaviside function. Its Legendre transform is 

/•I 2ti 
= 27r / dzP\z) = — — (P,+i(cos0^) - P;_i(cosr)) , (26) 

J cos 8" Zl + i 

the /-weighted scaled solid angle. 
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2. -Toy hat profile. We can skip the complications induced by the Legendre 
polynomials approximating the hat profile by 



with 



_ I «(Z+1) ic^^\ 

~\ fi'(-t)2 for > ^^^^ 



fi'^ = 27r(l-cose') . (2^ 



the scaled solid angle. 

3. -Gaussian profile. A large angle generalization of the well-known small 
angle Gaussian window: 

(Z+ "'" 

W^' = Q'e -2 — . (29) 



4.-SSSS profile. A profile of the form (e.g. ||T2 



W'(9,y)(x , (30) 



with some sort of cut off at 

1/2 



O'^iy) = arcsin ^min ^1, 



2y/^ - 1 



(31) 



the angle subtended by the horizon. We define this cut off with the Legendre 
transform: 

' iT^dH'e'^'^'l^'' for / < \\IQ^\ 
27rp£^^ for /> [1/0^ 

The hat profile is suggested by the belief that ^ circular averages are all 
that matters for 1 < I < 30. Circular averages have been extensively studied 
(eg. |]13|) and some useful properties have been found for them. The toy hat 
profile, due to its simplicity, is ideal for an heuristic explanation of the results 
obtained (Sec. |3.5|) . The Gaussian profile is particularly amenable to analytical 
work, thus providing a good check on the numerics. Furthermore it is probably 
a good approximation to the peaks found in P]. Finally, there is no harm in 
considering the SSSS profile, for it might have something to do with reality after 
all. 

With the W defined as they are, the functions C^{y) peak for different values 
of y. To facilitate comparison we should therefore define a reduced p^: 

Ps = = ^ = f (33) 

which ensures that, for high /, C^{y) peaks at the same ymax for the same for 
all profiles. For the SSSS profile C'(y) peaks always for y : d^{%j) ^ j. This is 
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Figure 2: C^^{y) for tlie liat profile (oscillating line), the toy hat profile (non- 
oscillating line), and Gaussian profile (with points). For all the curves = 0.2 
and < > is the same. 

understandable as the integrated anisotropy increases logarithmically as we go 
away from the core. As a result, the effective cross-section for SSSS textures is 
the horizon area, for all pg. For this reason we will always treat the SSSS profile 
separately, as it does not really have a parameter p^. 

It also happens that the value of C\ymax) depends on the profile, taking the 
form 

= (34) 

with a^"* ^ 0.58, a^^^ = 1, a^^'^^^ = l/(4e). By defining a reduced a as 

^ ^gauss 

a = av^^ = l.OSa'^'^* = 2a*°^ = (35) 
Ps 2v/^ 

we can ensure that C\ymax) for high I is the same for all the profiles with the 
same < > and p^. 

In Fig. we have plotted C^^iy) in units of < > n, with p^ = 0.2, for 
the first three profiles proposed. We notice that for y > ymax all the C\y) 
have the same form, as for all profiles W^^ cx (even though the proportionality 
constants may be different). Consequently, in this region C'(?/) oc fi*(?/) oc 1/N{y) 
regardless of the profile. For y <^ y-max we have W'"'^ ^ fi'^ and so C\y) is 
proportional to N{y) or smaller than that. The exact form of C\y) in this region 
is profile dependent, as it involves details of how W^^{fl^) starts to cut off at a 
particular scale. While this portion of the integrand is not the main contribution 
to the total C' it is clear that the detailed form of the spectrum is going to be 
profile dependent. It remains to be seen if these differences have any meaning 
once cosmic variance error bars are taken into account. 

It is also in the y -C ymax region that C'(?/) may be affected by the existence of 
non-spherical spot modes (e^ in (IH)). These are the only source of inaccuracies 
in the model proposed, and so whatever mistakes we make will be in this region. 
Like profile differences, spots' ellipticity can only affect the final C' in the fine 
detail. Most probably, ellipticity affects C' even less than the profile, as the spots 
seen in the simulations appear to be very circular. 

3.3 An expression for cr^(C;) 

We now turn to the issue of cosmic variance in the Ci. This turns out to depend 
on the nature of the theory. It is known that in Gaussian theories the Ci are 
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Figure 3: The relative non-Gaussian excess variance e'(ps) for Gaussian profile 
spots. 



X2/+i"distributed with a variance 

4W) = C"^^. (36) 

In some of the literature these are also the cosmic variance error bars attributed 
to non-Gaussian theories. We will use our model to compute directly (7^j^{Ci). 
From (|1^) we have that 

\ ' J mm' (n\n2n^n4) 

< Y^{n^,)Yl{n^,)Y^,{n^,)YU^^,) > . (37) 

The last factor can only be non-zero in one of the following cases: 

• if 77-1 = 712 and = but rii ^ n^, giving rise to [21 + 1)^ terms, 

• if 72i = 77-3 and n2 = (but rii ^ n2) or rii = and n2 = (but rii ^ 77,2), 
giving rise to a total of 2(2/ + 1) terms, 

• if ni = ri2 = ■'^s = "",4, giving rise to a single term. 
Putting all these terms together one has 

< >..^ ^ E < « > (1 27^) H- '-^1: < < > . 

(38) 

Subtracting off C'^ written as (^1]), using (p]), and neglecting terms in one 
finds after some algebra 

We note that the cosmic variance in the in texture scenarios is always larger 
than in Gaussian theories. We define e\ the relative size of the excess variance, 
as 

J _ ^'tx{Ci) - al{Q) ^^^^ 
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Figure 4: C spectra for the hat, toy hat, and Gaussian profiles. 



is a function of n, ps, and the profile. In general e\n,ps) = e'(l,ps)/n, so 
we define e'(ps) = e\l,ps), the excess variance per unit of 1/n. The quantity 
can be seen as a (theoretical) indicator of how non-Gaussian a given multipole 
is. A texture theory with n ^ oo will be a Gaussian theory. A low n theory 
will be very non-Gaussian, at least for some /. In Fig. | we have plotted e''{ps) 
for spots with a Gaussian profile. Generally the excess variance is very small for 
high / < 30. For low / (say from 2 to 5) the correction can be considerable, the 
more so the smaller the pg. 



3.4 C'*-spectra for different and profiles 

From (p^ we see that the parameters < > and n can only affect the overall 
normalization of the spectrum. Also ( |5^ suggests that the "average" spectrum 
normalization is proportional to irpl, although ps may affect the spectrum shape 
as well. Furthermore a roughly scale invariant (oc l/l^) spectrum can be expected. 
We do want more than a rough prediction, though, so we define a reduced spec- 
trum 

C- . J-'t^l'^^^ . (41) 



< >.2n\og{2).'7ip 



yl-k 



C describes in detail departures from scale invariance and factors out all that 
contributes only to the normalization of the spectrum. It depends only on the 
profile and p^. In Fig. ^ we have plotted the C'* spectrum for various values 
of Pg < 1/4 (to ensure causality), for the hat, toy hat, and Gaussian profiles. 
The spectrum normalization still depends on the profile. Even though the peak 
heights in C'(?/) have been matched by our definitions of p^ and a, the differences 
in the peak's shapes and in C'(?/) away from the peak produce different integrated 
for different profiles. Apart from this it is remarkable that the C'*(pg) spectra 
shapes do not depend on the profile in any significant way. For intermediate 
scales (/i < / < 30, for some /i), the spectra are scale-invariant if p^ is not too 
small. If Pg is small, the spectrum is slightly tilted. Whatever the case, for low 
/ (/ < /i) the power is suppressed relative to a scale invariant (or slightly tilted) 
spectrum. An exact evaluation of the tilt, /i, and low / suppression factor can 
only be done numerically. An approximate formula and heuristic explanation of 
this effect is given in Sec.|33. 
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Figure 5: C'* spectrum (middle line) and its Gaussian (lines) and non-Gaussian 
(bars) cosmic variance error bars for spots with a SSSS profile with n=0.25 and 
Ps = 0.2. 

Missing author supplied figures 

Figure 6: The spectral index n* for high / as a function of p^. 

3.4.1 SSSS collapses 

We have applied our formalism to concrete choices of {n,Pg) and profile which 
have often been suggested by simulations. In there is the suggestion that it is a 
good approximation to consider that only SSSS collapses followed by unwindings 
cause CMBR spots. This leads to a scenario with a low spot number density 
(n 0.25), and regardless of ps (which is a random variable in this case), to an 
effective spot size of the order of the horizon angular size. We have pointed out 
before that n and Ps do not matter for the spectrum shape produced by SSSS 
spots. We have plotted in Figure ^ the C'* spectrum {p^ = ps, a = a, no division 
by npl in (0)) for SSSS spots with ps = 0.2 (for definiteness) , with cosmic 
variance error bars computed from (|39|) . A flat spectrum without any significant 
low / cut-off is obtained, in agreement with This is a general feature for 

any profile whenever p^ is sufficiently large. The only novelty is the large non- 
Gaussian correction to the cosmic variance in this scenario, due to the low value 
of n (recent reruns of have shown an abnormally large cosmic variance). 

3.4.2 Gaussian spots 

In [§, on the other hand, it was found that both realistic (non-symmetric) un- 
windings and concentrations of gradient energy produce CMBR spots. Together, 
these two types of effect bring the spot number density up to n ^ 1. The 
spots produced by the two effects appear to be indistinguishable. Their profile is 
clearly non-SSSS, and is better approximated by a Gaussian profile with a small 
Pg, around 0.1. A clear evaluation of p^ has not been provided by the simula- 
tions. These simulations are essentially fiat space-time simulations (p) in which 
the horizon size is ambiguously defined. Hence we leave free parameter. 

Low Pg produce slightly tilted spectra even for intermediate / (cf. Fig. ^). A 
natural question is how the tilt depends on p^. We fitted the C'* spectra for 
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Figure 7: C'* spectrum (middle line) with its Gaussian (lines) and full (bars) 
cosmic variance error bars for Gaussian profile spots with = 0.25, 0.2, 0.1, and 
0.05, confronted with fitting Gaussian tilted spectra with the same normalization 
(points). 

Gaussian profiles to the infiationary type of spectra 



where is the spectral index. The fit was performed for I G (25,30), and the 
resulting n* = n^iPs) function is plotted in Fig. ^ For p^ > 0.15 we find ^ 1, 
but for pg = 0.05 and Pg = 0.1 we have respectively = 1.3 and n* = 1.2, 
for instance. Hence the importance of an accurate measurement for p^ as it will 
provide the intermediate scale spectral index for textures. Still, this is not the 
end of the story. In Fig. ^ we plotted the C'* spectrum for four values of p^ with 
1-0" cosmic variance error bars. Superposed on them are the fitting infiationary 
type of spectra, with their Gaussian cosmic variance error bars. If ps is not 
too small, a 1-a differentiation between the two theories arises for I = 2 and 3 
although the cosmic confusion at intermediate / is almost 1. This is somewhat 
surprising, as the cosmic variance in the low is very high (and even higher 
in texture scenarios). In spite of this, the extra suppression of low / in texture 
scenarios is strong enough to survive cosmic confusion. It is tempting to connect 
this effect with the abnormally low Ci for / = 2, 3 observed by COBE. However, 
the experimental error bars make it unwise to draw any conclusion. This effect 
also shows how a concept like the spectral index n*, coined for inflation, becomes 
inadequate as a texture spectra qualifler. Texture spectra have a non-uniform 
tilt. If Ps is very small the extra suppression, although meaningful if the cosmic 
variance were Gaussian (high n), is completely drowned by the non-Gaussian 
excess of cosmic variance. Although somewhat disappointing, this result also 
signals strong non-Gaussianity for low-/. At the same time as it renders the 
spectrum useless, it suggests that a particularly strong signal ought to exist in 
quantities measuring non-Gaussianity. 




(42) 
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3.5 An heuristic interpretation and an approximate for- 
mula 

The expansion into spherical harmonics acts as a scale filter. Each harmonic 
selects spots with an angular scale smaller than = jj^^j- Whatever the 

spot profile, if < then W^' ^ but if ft' > then W^' < ft'. 

In the latter case, depending on the profile, W''^ may either reach a plateau 
or simply decrease monotonically with il'. We can understand the 
shape of C''{y) by taking this into account. For y : Q,' < fi'*^*"*-*, each texture 
contributes to like Q*^, but when we add them up in the aj„ they may interfere 
either constructively or destructively. As a result, only a r.m.s. fluctuation 
proportional to -\/]V contributes to VC^. The overall contribution from a given 
y is therefore proportional to A^fi*^, which decreases as we go up in y. As y 
increases we do have more textures in the sky, but their apparent size is also 
much smaller, and since Nfl const, the balance makes C'(y) decrease. For 
y : Q' > the filter starts to act, collecting in a given a contribution from 

each individual texture proportional to VT'*^ ^ fi*^. This may either saturate or 
cut-off. Whatever the case the contribution from a given y will be proportional 
to N or less than that. The contribution from each texture is smaller than a 
constant, and C' increases as we go up in y, simply because there are more 
textures creating a r.m.s. fluctuation. Hence there will be a peak in C'(y) at 
y : ^'{y) ~ fi"'*"*^ and we may expect that most of the contribution to a given 
C' will come from the scales where the filter starts to act. The contribution from 
these textures is of the order A^f2^(*"*)^, and recalling that NQ' pa const, we can 
expect a value of C' proportional to oc This explains why we can 

always expect a roughly scale-invariant spectrum in the texture scenario. The 
only exception to this argument is the case where the Q' > regime does 

not exist (or is irrelevant), simply because is greater than < f2i >, the 

average angular size of the last texture. This is the case for low /, which do not 
act as a filter for any of the texture spots, because their cut-ofF scale is above the 
angular size of the last (and largest texture). Then, the largest contribution to 
C' comes from the last texture, this being the case for all I : ^< >. 

The last texture contribution to these C' is, of course, independent of /, and 
proportional to < fl^ >. Hence we can expect a white noise type of spectrum 
(C' oc const) in this I region. In practice, the transition from white noise to scale 
invariance happens very quickly. We never observe a white noise regime, but only 
a suppression for the low I. 

This heuristic argument can be converted into an approximate formula by 
taking the toy-hat profile and carrying out explicitly the integration, starting 
from y =< yi >. Let us define li as 
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Figure 8: Approximate spectra /(/ + 1)C' for various values of Ps- 
(the /-scale of the last texture), and ks as 

(the /-scale of textures living at the last-scattering surface). Then 

^ /-f^i W'^^ \ - ^is for / < /i 

^ ^ L "^^^ I Sfi'^^-^*) - 2^ - for / > h ^^^^ 

giving mathematical expression to what we have said above. This leads to the 
approximate spectral formula 

1 1 for / < /i 

for / > /i 



OC <; 3 _ 2/1^1+1) 1 f^^j^, l4bj 



1(1+1) + lUhs + 1) 

with /i ^ 0.5/ps and //« ~ for standard recombination, with ps = 2p^. 

The approximate spectra were plotted in Fig. 13 with an arbitrary vertical scale. 



Confronting this with Fig. |J, we see that (|46|) does provide a good qualitative 
description of the spectra. 

4 The two-point correlation function 

The two-point correlation function C{6) is a tool sadly closer to experiment than 
the C' spectrum (however see |T^ for a more hopeful point of view). We will 
now show that C (6) is a very bad discriminator between inflationary and texture 
scenarios. C{9) is defined as 

C{e) =< ^(fi)^(fi') >ots (47) 

with 6 = {Q,Q'), {Q and fl' any two fixed directions in the sky). Making use of 
(p!7D this can be written as 

c{e) = j:c^^p,{cose). (48) 

The theoretical C{9) can be estimated from the sky average of with 
a fixed 6 = {9^'): 

CoUO) = I ^^^-^in)^-^in')Sicos9 - cos ((]>)) (49) 
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which after some algebra leads to 

2/ + 1 

Cobs{d) = T.Ci^—Pi{cosO) (50) 

with Ci defined as in (|18|). An expression for C{9) in texture scenarios can be 
found by replacing (|2^ ) in (|48|) . It is curious that a simpler expression can be 
derived. Inverting the Legendre transform (^) one has 

9/ -I- 1 

W^{z,y)=Y.^W^\y)Pi{z). (51) 

If we define the Legendre-squared spot profile to be 

91 -I- 1 

W^^^\z, y) = J2^W^'\y)Piiz) (52) 

then the correlation function is simply 

C{9) = rdyNiy)W<'\cos9,y) . (53) 

47r Jo 

The profile W^^'^^ is the Legendre analogue of VT^T^riy^ (where -k is the convolution) 
in Fourier analysis. Depending on whether the profiles are defined in the 6 or 
in the Legendre spaces it may be easier to use ( [53| ) or to combine ( ^8]) and 
(|22D . If one wants to take into account the telescope beam filtering, one should 
multiply the C' in ( P5| ) by the square of the beam Legendre transform JF'^. In 
the numerics we will use an approximation to the COBE beam JF' = e'^C'+i^ 
with a = 27r(107360°). 

We have also derived an expression for the cosmic variance in Cobs{9)- In 
Gaussian theories (|^ with a factor of 2 correction) we have: 

^UCoUO)) = Y.^UQ) (^^j (cos^) (54) 

with aQ{Ci) given by (|36D . In the texture scenarios not only does one have to 
replace (Jq{Ci) by CTxi^i) (^s in (|39D), but also terms in PjP;/ (with / ^ I') 
appear, refiecting the inter-/ correlations. A rather elaborate calculation, in the 
moulds of Sec.pTB|, brings us to 

cr'MCobsm = E^''C'C''^^^^PKcos^)P,(cos0) (55) 

with 

25iv f a^(a^) \ J^'^dy N{y)W''^{y)W'''\y) 

"'~2/ + l^l^ ^ <a^ {Sy''dyN{y)W^'\y))Uy''dyN{y)W^''^{y)) ^ ' 
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Figure 9: t{9) for Gaussian spots witli {n = 1, = 0.25). 
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Figure 10: a^x (^op line), and its Gaussian (bottom line), diagonal (top points), 
and off-diagonal (bottom points) in Vw, contributions. 

where, again, we have neglected terms in 1/iVt. The intensity of the non-Gaussian 
correction to the cosmic variance can be measured by means of the quantity 

= ^m^) — 

which depends on n, p^, and the profile. The dependence on n is trivial, so 
we define e{ps,9) = e{n,ps, 9).n. The numerator in ( ^7|) can be rewritten as 
< (p^^(2))2 ^ showing that the variance in Cobs{9) is always larger in texture 
scenarios than in similar Gaussian theories. As before, we adopt the attitude 
that while this does reduce the predictability of the theory on C{9), it also signals 
non-Gaussian behaviour and the need of a non-Gaussian data-analysis approach 
to fully make out the predictions of the theory. In particular, on angular scales 
where e ^ 1, we know that it is worth studying the collapsed (to two points) 
n-point correlation function. In Fig. |^ we divided 9 G (0, vr) into 100 points 
and plotted e{9) for Gaussian spots with (n = 1, ps = 0.25). The fact that the 
low / multipoles contribute to C{9) for all 9 distributes non-Gaussianity nearly 
uniformly over 9. In Fig. |10] we analyse the various contributions to (t'^x^ sepa- 
rating its Gaussian, diagonal and off-diagonal components in Vw. We note that 
the inter-/ correlations responsible for the off-diagonal elements of Vui can act 
so as to reduce the cosmic variance, an important fact which we shall make use 
of in |Ty]. Naturally, the subtle differences in the C'-spectra pointed out at the 



end of Sec. |3.4| are completely drowned in C{9). As an example we have plotted 
in Fig. |TT] C {9) for the above texture theory, and for its fitting tilted spectrum 
inflationary theory. As C{9) spreads the low / over all 9, the two theories come 
out completely confused. Note however that the normalization in the C' and in 



C{9) can be substantially different (in Fig. |TT]we used a different normalization 
procedure than in Fig. 
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Figure 11: C{9) for Gaussian spots witli (n = 1, = 0.25) witli its full cosmic 
variance error bars, confronted with its fitting n = 1.001 theory (top line at 

5 The last texture and non-Gaussianity at low 

I 

In the derivation of (^) we made use of statistics for unordered textures. How- 
ever, a similar expression could have been obtained using ordered textures. In 
this case one should be wary for the W^'' are not independent variables. Never- 
theless formula (pT]) remains valid if one uses the marginal distribution functions 
of i/n in computing < PV^^^ > for each term in (0). The marginal distribution 
functions Pnijjn) are defined as 

PniVn) = J dyi...dyn-idyn+i.--dyNtPi...Ntiyi,--,yNt) (58) 
and using (p!0| ) they can be found to be 

Pniyn)=Niy^)e-'^(y-^^^^^ (59) 



where M{y) = Jq N{x)dx. Then: 

Att ^^Jo [n-iy. 



(60) 



and since Nt is very large we recover (^). The advantage of ( [501 ) is that it allows 
us to write 

= T.CL (61) 



with 

Cl = ^ [''dyr,N{y^)W^^\yr:)e-^^y-^^^^ (62) 

reporting the contribution of the n*^ texture to the multipole CK In Figure [T^ we 
plotted the C'* spectrum for four plausible scenarios: for a Gaussian profile with 
Ps = 0.05 and Ps = 0.1 (n = 1 in both cases), and for a SSSS profile with n = 0.25 
and n = 0.1 {ps = 0.2 for both). Underneath the spectra we plotted the result of 
stopping the sum ( |6TD at A^^ from 1 to 10, thus obtaining the contribution of the 
first 10 textures to each CK In all the cases considered we notice that the low 



18 



Missing author supplied figures Missing author supphed figures 
Missing author supphed figures Missing author supphed figures 



Figure 12: Contribution of the last ten textures for a Gaussian profile, with 
Ps = 0.05 and = 0.1 (n = 1 in both cases), and for a SSSS profile with 
n = 0.25 and n = 0.1 (and ps = 0.2). 

Missing author supplied figures Missing author supplied figures 
Missing author supplied figures Missing author supplied figures 

Figure 13: The and function for a Gaussian (top), and a SSSS profile 
(bottom). 



/ multipoles are dominated by the contribution from the last texture. The / at 
which this stops being the case depends on the spot details (profile, n, ps), but 
I = 2,3 always seem to be subject to this effect. Also the intensity of the last 
texture dominance at low I depends on the spot details, the effect being notably 
pronounced for Gaussian profiles with low ps- As we go up in /, the last texture 
becomes less prominent, but no other single texture replaces it as a dominant 
feature. We may ask how many textures are responsible for, say, 95% of a given 
for high /. Using the toy hat model and performing an approximate calculation 
similar to the one in Sec. |3.5| , we can show that this number is proportional to 
This implies that for high / the number of defects responsible for a given al^ is 
proportional to /. 

These results have far reaching consequences. If one believes in fiippant appli- 
cations of the central limit theorem, then one can expect the a^^ for high / < 30 
to be independent Gaussian distributed random variables, as they add up a large 
number of independent contributions (proportional to /). We have now shown 
that this argument can certainly not be repeated for the low /. In fact, our remark 
in Sec.^]3| that (TTxi^^i) > ^GaussiC^i)^ ^^e numerical results plotted in Fig. 
13, prove that the distinctly non-Gaussian at low /. This also implies that 

these aj„ are dependent random variables [|lT|. We can now understand better 
why this is so. The non- Gaussian features in the al^ for low / result from the 
fact that they are mostly due to one single lumpy object blatantly different from 
Gaussian noise. We may even expect the for low I to reproduce the morphol- 
ogy, size, and other features of the last defect. Devising a data- analysis method 
capable of uncovering these features is the purpose of [|lOl- We should stress that 
both the central limit theorem and Fig. ^ strongly suggest, but still do not prove, 
that the Gaussian for high / < 30. 

The non-Gaussian toy model studied in [|11| showed in which way the low / 
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multipoles would display non- Gaussian behaviour if they were fully due to the last 
texture. The contribution from n > 1 textures will soften this non-Gaussianity. 
Hence a plausible measure for the expected non-Gaussianity present in a low / 
multipole is a' = C{/C'', a quantity which depends on the profile, n, and ps. 
We have computed numerically and for a Gaussian and a SSSS profile. 
The results are plotted in Fig. |13|. A better connection between a' and detailed 
features of non-Gaussianity measure distribution functions will be given in []T0|. 



6 What's new? 

The focus in this paper was to derive for texture scenarios whatever one nor- 
mally derives for Gaussian theories (C', C{6), and their cosmic variances). We 
met anomalous behaviour (eg. formulae ( |39D and (|55D) which (theoretically) in- 
dicates non-Gaussianity, but we have left to the sequels and [0 the task 
of developing a non-Gaussian data analysis method. On the whole, within the 
standard lore, the results obtained are surprisingly similar to inflationary scenar- 
ios. Still, we found a few novelties. We discovered that texture spectra are tilted, 
but not uniformly. Fitting a tilted spectrum for / G (25, 30) (see Fig. P) leaves a 
significant suppression of power at low I G (2, 5) (see Fig. |^). We have also com- 
puted directly the cosmic variance error bars in Ci and Cobs{d)i and found them 



to be larger than in comparable Gaussian theories (see formulae (|39|) and (p5|) 
and Figs. |^ and |r^). In fact, if the non-Gaussian correction is very large (which 
does happen for most of the parameter values suggested in the literature), cosmic 
variance drowns the otherwise l-cr suppression of power at low / mentioned above 
(see Fig. However, one should not give way to despair. Large non-Gaussian 
corrections to cosmic variance only hint that we have not applied to the theory 
the right data analysis procedure. To give a flavour of [1^ let us point out that 



as the Ci are dependent random variables in non-Gaussian theories ([|TT]), it may 
happen that the Ci spectrum shape seen by any observer is never the average C' 
spectrum shape. A more intelligent method to make predictions on Ci spectra 
and to estimate global parameters in the C' spectrum is in order. It appears that 
even a proper study of the C' (or Ci) spectrum for textures requires non-Gaussian 
data analysis. 

The central result in this paper is undoubtedly the proof that texture low / 
multipoles are strongly non-Gaussian, and that the last texture is to be blamed 
for this (see Sec.^. We showed how the low C' are mostly due to the last texture 
(slightly perturbed by the one after - see Fig. |l^), and so we may expect the 
last texture's non-Gaussian features to be imparted on the low / multipoles. In 



11] we proposed the use of m-structure and inter-/ measures to complement 
the CK m-structure measures act as multipole shape factors. Inter-/ measures 
correlate preferred directions in two multipoles. These two types of quantities 
are uniformly distributed in Gaussian theories. In topological defect scenarios we 
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can expect low / shape factors to reflect the morphology of the last defect. Hence 
their distributions should be non uniform, peaking at different values for, say, 
texture and cosmic strings. Inter-/ measures should also display the correlations 
between the various low I which are all due to the same last defect. In remains 
to be seen ([jlO[) how much cosmic confusion there is between these signals and 
their Gaussian counterparts. 

On a lower key this paper was simply intended as an analytical model for 
texture CMBR skies with the greatest possible generality. We tried to leave as 
free parameters whatever simulations have failed to determine. The point is to 
find out what exactly must be decided in order to answer a particular question. 
We have concluded, for instance, that the C' spectrum shape is insensitive to n or 
a, and depends very little on the spots' profile, once the identifications (^) and 
(|35D have been made. The controlling parameter is which fixes the spectrum 
tilt at intermediate / (see Fig. |^). The spectrum normalization, on the other 
hand, depends on all the parameters of the model (see Sec. p.l|) , which must be 
evaluated before any prediction for the texture symmetry breaking energy or 
bias parameter can be made. Interestingly enough, we found that the low I non- 
Gaussianity effects depend only on n, p^, and profile (Sec.|^). These are not free 
parameters of the theory, so one ought to decide on this important issue purely 
by means of hard work. 

The paradigm proposed here cannot easily be adapted to cosmic string sce- 
narios. Neglecting non-circular spot modes is, in texture scenarios, only slightly 
inaccurate. The string counterpart to circular spots, line discontinuities cut-off at 



horizon distance (as in [|T6[), is definitely a much cruder approximation. More im- 
portantly, neglecting texture spot-spot correlations introduces only a small error 
(in fact it would not be difficult to incorporate small correlations in our model). 
On the contrary, the various segments of the string network form a Brownian 
random walk. Segments linked to each other are therefore very highly correlated, 
making an analytical treatment cumbersome. Uncorrelated string segments (as 
in are not very realistic. Overall, although we can boast a more realistic set 



of assumptions than |]T6| , this is mostly due to the fact that textures are far sim- 
pler than strings. This also justifies why we can go further than |T^. We treated 
with greater rigour (no "ergodic" assumptions) a much broader set of models. 
We also managed to treat exactly (within the model) the problem of cosmic vari- 
ance. Finally, greater sophistication in the statistical machinery opened doors to 
the issue of non-Gaussian data analysis in defect scenarios (developed further in 

mm). 
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Abstract 

We propose an analytic method for predicting the large angle CMBR 
temperature fluctuations induced by model textures. The model makes use 
of only a small number of phenomenological parameters which ought to be 
measured from simple simulations. We derive semi- analytically the C'- 
spectrum for 2 < / < 30 together with its associated non-Gaussian cosmic 
variance error bars. A slightly tilted spectrum with an extra suppression 
at low / is found, and we investigate the dependence of the tilt on the 
parameters of the model. We also produce a prediction for the two point 
correlation function. We flnd a high level of cosmic confusion between tex- 
ture scenarios and standard inflationary theories in any of these quantities. 
However, we discover that a distinctive non-Gaussian signal ought to be 
expected at low /, reflecting the prominent effect of the last texture in these 
multipoles. 
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1 Introduction 



In the recent past, a large amount of work has been directed towards predicting 
the CMBR temperature anisotropics associated with inflationary scenarios (see 
[1] for a review), and with the various topological defect scenarios (e.g. [2, 3]). 
Inflationary scenarios are far better explored in this respect. The main reason 
for this is that it is easier in these models to plead ignorance of the detailed 
mechanism responsible for the density fluctuations of the Universe. Quantum 
fluctuations in the metric are a fuzzy subject, and hand-waving arguments for 
Gaussian fluctuations with a particular type of spectrum are normally considered 
acceptable. Once this is done, one is free to take advantage of the well developed 
and sophisticated industry of methods aimed at treating Gaussian fluctuations 
[4]. A precise prediction for the outcome of concrete experiments follows easily, 
making inflation popular with experimentalists. The state of affairs in topological 
defect scenarios is rather different. No lack of fundamental physics understanding 
hampers working out in rigour the defect network evolution and the fluctuations 
they induce on the matter and radiation of the Universe. The price to pay for 
our honesty is that a serious treatment of the problem is computationally over- 
powering. A reliable link between the defect network and the perturbations it 
induces is undoubtedly still missing. To top off the trouble, it turns out that 
the fluctuations induced by defects are often non-Gaussian. Setting up a data- 
analysis framework geared towards non- Gaussian fluctuations is virgin ground 
(but see [4, 5, 6, 7] for Gaussianity tests). Overall, even though defect scen- 
arios have so far been poor in solid predictions, they have a potentially polemic 
value not to be discarded. Non-Gaussian data analysis does not exist partly be- 
cause non-Gaussian theories are too vast a class to be treated in one go. Defect 
scenarios provide a motivated and topical example with which one can rehearse 
non-Gaussian data analysis techniques and experimental strategies. 

The main technical complication in predicting CMBR fluctuations in defect 
scenarios arises from the fact that ^ maps photograph the whole sky (that is, 
our past light cone clipped by the last scattering surface) and not just its inter- 
section with the last scattering surface. Hence the defect network, surrounding 
matter, and CMBR photons have to be evolved through a large number of ex- 
pansion times. The problem is usually tackled by means of full-blown simulations 
[2, 3]. These require access to supercomputers, and even so face severe limita- 
tions. Important issues like cosmic variance in defect scenarios have hardly been 
addressed. A possible computational short cut can be found if one is to assume 
scaling. The defects network is expected to scale, by which one means that it 
looks statistically the same at any time, once the network length scale is equated 
to the horizon distance. This suggests dividing the sky cone into cells correspond- 
ing to expansion times and horizon volumes. Scaling makes it plausible that the 
CMBR pattern in each of these cells is statistically the same, once some angular 
scaling procedure is applied to it. This being the case, all that is required from 
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the simulations is good statistics for the CMBR spots induced by defects inside 
horizon size boxes during one expansion time. Such simulations are considerably 
easier to perform. A prototype of this type of study is [8]. In the case of tex- 
tures, the short-cut which we have described has been adopted by [9] where a 
"scaling-spot-throwing" process was implemented on a computer. The spots were 
all derived from model SSSS (self-similar and spherically symmetric) collapses. In 
this paper we advocate the use of analytical techniques for scaling-spot-throwing. 
An analytical approach allows the derivation of exact formulae for the C\ C{9), 
and their cosmic variance, from which interesting general properties can be de- 
rived. These properties are naturally less evident in computer simulations, unless 
one explicitly sets out to find them. Furthermore, an analytical approach allows 
the treatment of a more general case than [9]. The method proposed here applies 
to quasi-circular spots with any profile. The scaling size (ps), number density (n), 
and brightness factor (a) of the spots are also left as free parameters. We should 
admit that only the brightness factor is a truly free parameter in the texture 
model. All others should be considered as phenomenological in origin, eventually 
to be fixed by simulations. As they still constitute a controversial matter, we 
choose to consider a broader class of model spots. We can then sort out which 
spot properties it is really important to measure from the simulations. 

The plan of this paper is as follows. In Sec. 2 we start by deriving the statistical 
properties of scaling texture model collapses as they appear in the sky cone, and 
from these a formal solution for the joint distribution function of the CMBR 
is given. Using this solution, in Sees. 3.1 and 3.3 we derive exact formulae for 
the C' and c(C;) in texture models. We prove the general result that the cosmic 
variance in the C; is larger in texture scenarios than in Gaussian theories. The 
excess variance is negligible for high /, but significant at low /. In Sees. 3. 2 and 
3.4 we compute numerically the C' spectra (with cosmic variance error bars) for 
a large class of spot types. We find that the spots' profile, intensity, and number 
density affect only the spectrum normalization. The spectrum shape is largely 
controlled by the spots' scaling size ps alone. In general, the spectra are slightly 
tilted (and we numerically compute the dependence of the tilt on p^), but with 
an extra suppression of power at low /. In Sec. 4 we derive expressions for the 
two-point correlation function C{6) and its cosmic variance. Although closer to 
experiment, C{6) appears to be a very bad discriminator between texture and 
infiationary scenarios. The most original result proved in this paper is presented 
in Sec. 5, where we ask how many textures are responsible for a given a^. For 
high /, this number goes like /, suggesting Gaussian behaviour. For low / we find 
that the are mostly due to the effects of a single texture: the last (that is, 
the closest) texture. This ought to impart a peculiar non-Gaussian signal to the 
low a^, a fact which we prove explicitly. The purpose of the sequel to this paper 
[10] is to devise a low-/ data-analysis technique (in the moulds of [11]) capable of 
fully characterising this effect. 
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2 Analytical spot-throwing 



In topological defect scenarios the CMBR fluctuations at large angles are due 
only to the ISW (integrated Sachs- Wolfe) effect caused by the defect network 
([2]). For a standard recombination scenario the validity range for this approx- 
imation is 2 < / < 30. If early reionization occurs, the upper limit for this range 
becomes the value of / corresponding to the angular scale of the horizon at "last 
scattering" (defined in some conservative way). As it is not clear which is the 
most sensible scenario we take as a working hypothesis the choice which maxim- 
izes the applicability of our model: standard recombination. Now if we assume 
that the statistical scaling of the texture fields extends to their associated time- 
dependent metric perturbations, then we may expect the ISW spots themselves 
to scale. This allows the direct modelling of the statistics for texture induced 
spots as they appear in the sky. From these, we can then write an expression for 
the joint distribution function for 2 < / < 30. For simplicity, we will assume 
that the texture spots are independent (not required by the formalism, but a 
good approximation). Then each texture angular position in the sky is uniformly 
distributed. Furthermore its time position in the sky cone can be associated with 
a modified Poisson process in time, in which the probability of an event is allowed 
to vary. Let n be the average number of texture spot producing configurations 
per horizon volume. Then their average volume density is 

m 



97rt3 



and using the relation between the proper area dS^ the time of emission t, and 
the apparent solid angle JO: 





ds = dn yst I 

we find that the angular density of spots in the sky per unit of emission time is 

(3) 

It is convenient to introduce the variable y = log2(to/t), the expansion time 
number as we go backwards in time. Fach Ay = 1 represents one expansion 
time, starting from j/o = (here and now) and going as far as yis = log2(to/t;s), 
when we hit the last scattering surface. In the standard recombination scenario 
yig 15. In terms of y the spot density in the sky is 

N{y) = 2nlog{2){2y/'- if . (4) 
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Let us now introduce a Poisson-type process in j/, with a variable probability 
density N[y). If pn{y) is the probability of having exactly n textures at y' < j/, 
then 

^-E^ = Niy)ip^_,-p^) , p_,=0. (5) 

This hierarchy of equations fully specifies the statistics of the textures' time- 
position in the sky. From the Pn{y) one can also derive the probability density of 
the n^^ texture position j/„ by means of 

Pn{yn) = Pn-l{yn)N{yn) . (6) 

However one should be wary of the j/„ variables, as they are not independent. Not 
only are their ranges constrained by j/i < j/2 < ••• < J/n but also their statistical 
dependence is patent in the conditional distribution functions. From (5) and (6) 
one has that 

P,{y,) = N{y^)e-fo''dy^(y) (7) 

but then 

P2(j/2|j/i) = iV(y,)e-^"^^^(^) (8) 

which is dependent on j/i. From these expressions we can derive the joint distri- 
bution function 

^12(2/1,2/2) = iV(j/i)iV(j/2)e-^o"^^M^) (9) 
which misleadingly factorizes. More generally one has 

Pi....(?/i,?/2,...,J/.) = iV(j/i)iV(j/2)...iV(j/.)e-r"^^^(^) . (10) 

The statistics contained in P„ or pn refer to a sequence of textures ordered in y. 
Ordering will be crucial in Section 5 when we uncover the last texture. However 
for some other purposes (like the derivations in Section 3) we may simply consider 
unordered textures. Then we have a total of Nf = Jq'' dy N[y) ^ 1 textures in 
the sky. Their positions j/„ are independent random variables with a distribution 

PniVn) = ^ (11) 

for all n, which is considerably simpler. 

We now model the texture spot patterns. These ought to follow some sort of 
statistical scaling law, where the scale angle is 

e'{y) = arcsin (^min (^1, (12) 

for a pattern laid at time y. Here ps is the scaled impact parameter [ps = pj {'it)) 
giving (Tg = 7rp^, the scaled cross section for photon anisotropics. Causality 
requires that p^ < 1/2. It may happen that p^ is a random variable itself, but 
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we shall ignore this complication. For small, nearly circular spots with profile 
W^{6^y)^ scaling implies that 

W'{e,y) =W{x) (13) 

with X = 6/6^. This is a sensible scaling law even for j/ < 1, although for y < I 
the exact scaling law should be more complicated, as the photons propagating 
through the scaling metric no longer move along parallel trajectories. In partic- 
ular, if J/ < j/o = 31og(l + Ps), we live inside a texture and are enclosed by its 
pattern. Unless {2nps)^ 1 textures like these are very rare, and the approxim- 
ation (13) should produce good enough final results. 

We seek to write the coefficients for skies filled with nearly circular spots 
described by the above statistics. For each of these spots let us first point the 
2;-axis at their centre, and perform the spherical harmonic decomposition in that 
frame. We obtain an expression of the form 



hlin) = a.^W^'{y.^).^^.6.^o + 4 (14) 

where W'"^' is the Legendre transform of the scaled profile: 

= 27r y"'^ dz P'{z)W'{z, yn) (15) 

and is a perturbation induced by the eccentricity. This perturbation is negli- 
gible, even if the spots are not very nearly circular, but / <^ [l/^*(j/„)]. The extra 
factor ttn in (15) accounts for the texture brightness, and is a random variable 
with a j/„-independent distribution function. Performing a rotation to a general 
frame, where the n^^ texture has coordinates 0„, and summing over all textures 
one finally obtains: 

aL = j:<^nW^'iyn)Y::{n^). (16) 

n 

The variables a„ are independent but equally distributed. Their distribution has 
to be determined from the simulations. The 0„ are uniformly distributed, and 
from them one can determine the 1^(0^) distributions (but with care, as the Yj^ 
are not independent). Finally from (10), (12), and (15) one can determine the 
W^' distributions. The W^' makes the various terms in (16) dependent variables, 
rendering an analytic approach for the distributions through (16) unfeasible. 
In any case finding the marginal distributions of the would not be the end 
of the story, since the are necessarily dependent random variables (see [11]). 
Hence we still would have to find the joint distribution F(a^^, a^^' •••) ^ com- 
plete solution to the problem of cosmic variance in texture scenarios. Fxpression 
(16) constitutes a formal solution to this problem which, despite all the obstacles 
to an analytical solution, is particularly suited to Monte Carlo simulations. It 
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separates the individual effect of eacli texture on eacli of tlie a^. It also factor- 
izes each texture contribution into three factors. The are purely geometrical 
factors which should be present in any SSS (statistically spherically symmetric) 
theory (see [11] for a definition). The W^' are a refiection of the structure of the 
sky (foliation into expansion times and horizon volumes) and should be present 
whenever there is scaling. Finally, the a„'s factor out all that is peculiar to the 
texture model under consideration, and will have to be measured from detailed 
texture simulations. 



3 The angular power spectrum 
3.1 An expression for the 

Throughout this paper we will use the notation C; for the angular power spectrum 
of realizations, and C' for its ensemble average. Whereas C; is a random variable, 
C' is a number. In any SSS theory one has 

< ala^::, C'8\,8Z, . (17) 

The C' can be estimated from the observable 

^ — _/ 



since < C; >obs= C' in any SSS theory. Cosmic variance introduces error bars of 
size cr[Ci) into this estimate. We start by deriving an exact expression for the 
C'. From (16) and (18) we have that 

C =< Ci >= < ^nan' X wlwl X yj:(o.)yj,(o.o > (19) 

™ nn' 

but since 

< Yi*{n^)Yi{n^,) ^ (20) 

47r 



we have simply 



C = <Wi^> . (21) 

47r 



Computing the averages using the unordered statistical description (11) we can 
write 

C' = —r^ dyN{y)W'^'{y). (22) 
47r Jo 

lips is a random variable with distribution f{ps)^ then expression (22) still holds 
but now with W'"^'"^ replaced by / dps f {ps)W''^'^ (y , Ps) ■ Fxpression (22) lends itself 
to be written as 

C'= / dyC'iy) (23) 
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Figure 1: The density C\y) multiplied by /(/ + 1) for / from 2 to 30 for (n = 
= 0.4) and (n = 1,^^ = 0.08) textures with Gaussian profile spots. 

thus defining a density 



measuring the contribution of textures living at time j/ to a given C . In Fig. 1 
we have plotted /(/ + l)C\y) for (n = 1, = 0.4) and (n = 1, = 0.08) textures 
with Gaussian profile spots, for / from 2 to 30. These plots are to be confronted 
with the common belief that a given C' is mostly due to textures living at a time 
when the apparent size of the horizon equals the multipole angular scale. We 
find instead that C\y) peaks for y: 0*(j/) j?^^, that is, when the textures' 
apparent size fits the multipole angular scale. Furthermore the peaks comprise 
several generations and an accurate prediction for the C' spectrum can never 
neglect any of the integrand in (23). Changing n only affects the normalization 
of C\y). Reducing not only decreases the normalization but also shifts the 
peaks to the left necessarily affecting the shape of the integrated CK In general 
as / increases the peak's heights (multiplied by /(/ + I)) are slightly reduced, but 
their width increases. The delicate balance of these two effects will determine 
whether the C' spectrum is fiat or not. 

3.2 Profiles and reduced a and ps 

It is not clear what profile W real texture spots have. Simulations [2, 8] suggest 
that the well-known SSSS profile may be oversimplified. We find that the exact 
profile form has little impact on the C'-spectra. For intermediate /, profiles do 
not matter at all, and for low / they induce differences easily confused by cosmic 
variance. The issue of matching free parameters for different profiles is, however, 
non-trivial. We have considered the following profiles: 




(24) 
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1. -Hat profile. Defined by 

W\x) = H{1- x) (25) 
where H is tlie Heaviside function. Its Legendre transform is 

W'' = 27r C dz P'{z) = -^{Pi+^{cos9') - Pi_^{cos9'j) , (26) 
J cose-' 2/ + i 

tlie /-weiglited scaled solid angle. 

2. -Toy hat profile. We can skip the complications induced by the Legendre 
polynomials approximating the hat profile by 

,„ f for < ^l{sat) ^ 4n 

\ 0'(-t)2 for 0« > O'(-t) ^^'^ 

with 

0^ = 27r(l -cos^) . (28) 

the scaled solid angle. 

3. -Gaussian profile. A large angle generalization of the well-known small 
angle Gaussian window: 

W'' = n'e . (29) 

4. -SSSS profile. A profile of the form (e.g. [12]) 

W'(e,y)oc , ^ (30) 



wi 



ith some sort of cut off at 



9''{y) = arcsin (^mm (^1, ^^ji^j j (31) 

the angle subtended by the horizon. We define this cut off with the Legendre 
transform: 

' 27r9''9'e-'^^'/^' for / < [1/9''^ 
2Tr9'^Z^ for/>[l/^^ 



^^'= ^^^V (32) 



The hat profile is suggested by the belief that ^ circular averages are all 
that matters for 2 < / < 30. Circular averages have been extensively studied 
(eg. [13]) and some useful properties have been found for them. The toy hat 
profile, due to its simplicity, is ideal for an heuristic explanation of the results 
obtained (Sec. 3. 5). The Gaussian profile is particularly amenable to analytical 
work, thus providing a good check on the numerics. Furthermore it is probably 
a good approximation to the peaks found in [8]. Finally, there is no harm in 
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considering the SSSS profile, for it miglit liave sometliing to do witli reality after 
all. _ 

With the W defined as they are, the functions C\y) peak for different values 
oi y. To facilitate comparison we should therefore define a reduced p^: 

hat toy 

p=^ = ^= pa^^'' (33) 
1.8 2 ^ ^ ' 

which ensures that, for high /, C\y) peaks at the same ymax for the same p^ for 
all profiles. For the SSSS profile C\y) peaks always for y : 0^{y) y. This is 
understandable as the integrated anisotropy increases logarithmically as we go 
away from the core. As a result, the effective cross-section for SSSS textures is 
the horizon area, for all p^. For this reason we will always treat the SSSS profile 
separately, as it does not really have a parameter p^. 

It also happens that the value of C'{ymax) depends on the profile, taking the 
form 

rh.. \ < > •2nlog(2).7rpg ^ 

C (ymax) = 1(1+1) ^ ^ 

with a'"'* 0.58, = 1, a^""** = l/(4e). By defining a reduced a as 

gauss 

a = ay^^ = 1.05a''"* = 2a*°^ = ^ (35) 

Ps " ' 



we can ensure that C'(ymax) for high / is the same for all the profiles with the 
same <a^> and p^. 

In Fig. 2 we have plotted C^^(y) in units of < > n, with p^ = 0.2, for 
the first three profiles proposed. We notice that for y > ymax all the C\y) 
have the same form, as for all profiles W^^ oc 0* (even though the proportionality 
constants maybe different). Consequently, in this region C\y) oc 0*(j/) oc 1/N(y) 
regardless of the profile. For y <^ ymax we have W''"^ <^ 0* and so C\y) is 
proportional to N(y) or smaller than that. The exact form of C\y) in this region 
is profile dependent, as it involves details of how W^^(^^) starts to cut off at a 
particular scale. While this portion of the integrand is not the main contribution 
to the total C' it is clear that the detailed form of the spectrum is going to be 
profile dependent. It remains to be seen if these differences have any meaning 
once cosmic variance error bars are taken into account. 

It is also in the y <^ ymax region that C\y ) may be affected by the existence of 
non-spherical spot modes (e^ in (14)). These are the only source of inaccuracies 
in the model proposed, and so whatever mistakes we make will be in this region. 
Like profile differences, spots' ellipticity can only affect the final C' in the fine 
detail. Most probably, ellipticity affects C' even less than the profile, as the spots 
seen in the simulations appear to be very circular. 
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Figure 2: C^^{y) for the hat profile (oscillating line), the toy hat profile (non- 
oscillating line), and Gaussian profile (with points). For all the curves = 0.2 
and < > is the same. 



3.3 An expression for <t^(C/) 

We now turn to the issue of cosmic variance in the C;. This turns out to depend 
on the nature of the theory. It is known that in Gaussian theories the C; are 
X2/+i"di^t^it)uted with a variance 

= C"^^ . (36) 

In some of the literature these are also the cosmic variance error bars attributed 
to non-Gaussian theories. We will use our model to compute directly (Tj-^{Ci). 
From (16) we have that 

<Cf>obs = (21+lY ^ ^ < «m «-^2 «-^3 «-^4 >< > 

^ ^ mm' (nin2n3rn) 

<y:*(o.jy:(o.jy:i(o.3)y:,(o.j> . (37) 

The last factor can only be non-zero in one of the following cases: 

• if ni = ^2 and = ^4 but rii ^ ns, giving rise to (2/ -|- 1)^ terms, 

• if ni = ns and ^2 = ^4 (but rii ^ ^2) or rii = ^4 and ^2 = (but rii ^ ^2), 
giving rise to a total of 2(2/ -|- 1) terms, 

• if ni = ^2 = ^3 = ^4, giving rise to a single term. 
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Figure 3: The relative non- Gaussian excess variance e'{ps) for Gaussian profile 
spots. 



Putting all these terms together one has 



(38) 

Subtracting off C'^ written as (2f), using (f f), and neglecting terms in 1/Nt, one 
finds after some algebra 

We note that the cosmic variance in the C' in texture scenarios is always larger 
than in Gaussian theories. We define e', the relative size of the excess variance, 
as 

e' is a function of n, p^, and the profile. In general e'{n,ps) = e'(l,ps)/n, so 
we define e'{ps) = e'(l,ps), the excess variance per unit of 1/n. The quantity e' 
can be seen as a (theoretical) indicator of how non-Gaussian a given multipole 
is. A texture theory with n — )• oo will be a Gaussian theory. A low n theory 
will be very non-Gaussian, at least for some /. In Fig. 3 we have plotted e'(ps) 
for spots with a Gaussian profile. Generally the excess variance is very small for 
high / < 30. For low / (say from 2 to 5) the correction can be considerable, the 
more so the smaller the ps. 
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Figure 4: C * spectra for the hat, toy hat, and Gaussian profiles. 

3.4 C'*-spectra for different Pg and profiles 

From (24) we see that the parameters < > and n can only affect the overall 
normalization of the spectrum. Also (34) suggests that the "average" spectrum 
normalization is proportional to it pi, although ps may affect the spectrum shape 
as well. Furthermore a roughly scale invariant (oc l/P) spectrum can be expec- 
ted. We do want more than a rough prediction, though, so we define a reduced 
spectrum 

&^ = ^"^(^+^)^' (41) 
< a2 >.2nlog(2).7rp2 ' ^ ^ 

C'* describes in detail departures from scale invariance and factors out all that 
contributes only to the normalization of the spectrum. It depends only on the 
profile and p^. In Fig. 4 we have plotted the C'* spectrum for various values 
of p^ < 1/4 (to ensure causality), for the hat, toy hat, and Gaussian profiles. 
The spectrum normalization still depends on the profile. Fven though the peak 
heights in C\y ) have been matched by our definitions of p^ and a, the differences 
in the peak's shapes and in C\y ) away from the peak produce different integrated 
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Figure 5: C'* spectrum (middle line) and its Gaussian (lines) and non-Gaussian 
(bars) cosmic variance error bars for spots with a SSSS profile with n=0.25 and 
Ps = 0.2. 

C' for different profiles. Apart from this it is remarkable that the C'^iPg) spectra 
shapes do not depend on the profile in any significant way. For intermediate 
scales (/i < / < 30, for some /i), the spectra are scale-invariant if is not too 
small. If p^ is small, the spectrum is slightly tilted. Whatever the case, for low 
/ (/ < /i) the power is suppressed relative to a scale invariant (or slightly tilted) 
spectrum. An exact evaluation of the tilt, /i, and low / suppression factor can 
only be done numerically. An approximate formula and heuristic explanation of 
this effect is given in Sec. 3. 5. 

3.4.1 SSSS collapses 

We have applied our formalism to concrete choices of (n^p^) and profile which 
have often been suggested by simulations. In [9] there is the suggestion that it is a 
good approximation to consider that only SSSS collapses followed by unwindings 
cause CMBR spots. This leads to a scenario with a low spot number density 
(n 0.25), and regardless of ps (which is a random variable in this case), to an 
effective spot size of the order of the horizon angular size. We have pointed out 
before that n and ps do not matter for the spectrum shape produced by SSSS 
spots. We have plotted in Figure 5 the C'* spectrum (p^ = p^, a = a, no division 
by Trpl in (41)) for SSSS spots with ps = 0.2 (for definiteness), with cosmic 
variance error bars computed from (39). A fiat spectrum without any significant 
low / cut-off is obtained, in agreement with [2, 9]. This is a general feature for 
any profile whenever p^ is sufficiently large. The only novelty is the large non- 
Gaussian correction to the cosmic variance in this scenario, due to the low value 
of n (recent reruns of [2] have shown an abnormally large cosmic variance). 
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Figure 6: The spectral index lor high / as a lunction ol p^. 
3.4.2 Gaussian spots 

In [8], on the other hand, it was found that both realistic (non-symmetric) un- 
windings and concentrations ol gradient energy produce CMBR spots. Together, 
these two types ol effect bring the spot number density up to n 1. The 
spots produced by the two effects appear to be indistinguishable. Their profile is 
clearly non-SSSS, and is better approximated by a Gaussian profile with a small 
around 0.1. A clear evaluation ol has not been provided by the simula- 
tions. These simulations are essentially fiat space-time simulations ([8]) in which 
the horizon size is ambiguously defined. Hence we leave p^ as a tree parameter. 
Low p^ produce slightly tilted spectra even for intermediate / (cl. Fig. 4). A 
natural question is how the tilt depends on p^. We fitted the C'* spectra for 
Gaussian profiles to the infiationary type ol spectra 

n' —l \ I 9 — n' 



r / + ^ r 



C' oc -) ' .( ) ' .( (42) 



where is the spectral index. The fit was performed for / G (25,30), and the 
resulting = n^{Pg) lunction is plotted in Fig. 6. For p^ > 0.15 we find 1, 
but for p^ = 0.05 and p^ = 0.1 we have respectively = 1.3 and = 1.2, 
for instance. Hence the importance ol an accurate measurement for p^ as it will 
provide the intermediate scale spectral index for textures. Still, this is not the 
end ol the story. In Fig. 7 we plotted the C'* spectrum for four values ol p^ with 
l-cr cosmic variance error bars. Superposed on them are the fitting infiationary 
type ol spectra, with their Gaussian cosmic variance error bars. II ps is not 
too small, a l-cr differentiation between the two theories arises for 1 = 2 and 3 
although the cosmic contusion at intermediate / is almost 1. This is somewhat 
surprising, as the cosmic variance in the low C' is very high (and even higher 
in texture scenarios). In spite ol this, the extra suppression ol low / in texture 
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ps=0.25 and fitting n=1.001 ps=0.2 and fitting n=1.06 
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Figure 7: C'* spectrum (middle line) with its Gaussian (lines) and full (bars) 
cosmic variance error bars for Gaussian profile spots with = 0.25, 0.2, 0.1, and 
0.05, confronted with fitting Gaussian tilted spectra with the same normalization 
(points). 

scenarios is strong enough to survive cosmic confusion. It is tempting to connect 
this effect with the abnormally low C; for / = 2, 3 observed by COBE. However, 
the experimental error bars make it unwise to draw any conclusion. This effect 
also shows how a concept like the spectral index n\ coined for infiation, becomes 
inadequate as a texture spectra qualifier. Texture spectra have a non-uniform 
tilt. If ps is very small the extra suppression, although meaningful if the cosmic 
variance were Gaussian (high n), is completely drowned by the non-Gaussian 
excess of cosmic variance. Although somewhat disappointing, this result also 
signals strong non-Gaussianity for low-/. At the same time as it renders the C' 
spectrum useless, it suggests that a particularly strong signal ought to exist in 
quantities measuring non-Gaussianity. 
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3.5 An heuristic interpretation and an approximate for- 
mula 

The expansion into spherical harmonics acts as a scale filter. Each harmonic 
selects spots with an angular scale 0* smaller than 0'(««*) = Whatever the 

spot profile, if 0* < 0'(*"*) then W'' ?^ 0% but if 0* > 0'(*"*) then W'' < 0^ 
In the latter case, depending on the profile, W''"^ may either reach a plateau 
(sa 0'(*"*)) or simply decrease monotonically with 0*. We can understand the 
shape of C\y) by taking this into account. For j/ : 0* < 0'(*"*), each texture 
contributes to C' like 0*^, but when we add them up in the they may interfere 
either constructively or destructively. As a result, only a r.m.s. fiuctuation 
proportional to \/iV contributes to a/ZT'. The overall contribution from a given 
y is therefore proportional to A^O*^, which decreases as we go up in j/. As j/ 
increases we do have more textures in the sky, but their apparent size is also 
much smaller, and since A^O const, the balance makes C\y) decrease. For 
J/ : 0* > 0'(*"*), the filter starts to act, collecting in a given C' a contribution from 
each individual texture proportional to W^^'^ <^ 0*^. This may either saturate or 
cut-off. Whatever the case the contribution from a given y will be proportional 
to A^ or less than that. The contribution from each texture is smaller than a 
constant, and C' increases as we go up in y, simply because there are more 
textures creating a r.m.s. fiuctuation. Hence there will be a peak in C\y) at 
y : 0*(j/) 0'(*»*), and we may expect that most of the contribution to a given 
C' will come from the scales where the filter starts to act. The contribution from 
these textures is of the order A^0'(*"*)^, and recalling that A^O* sa const, we can 
expect a value of C' proportional to 0'(««*) qc l/P. This explains why we can 
always expect a roughly scale-invariant spectrum in the texture scenario. The 
only exception to this argument is the case where the 0* > 0'(««*) regime does 
not exist (or is irrelevant), simply because 0'(««*) greater than < Oi >, the 
average angular size of the last texture. This is the case for low /, which do not 
act as a filter for any of the texture spots, because their cut-off scale is above the 
angular size of the last (and largest texture). Then, the largest contribution to 
C' comes from the last texture, this being the case for all / : 0'(««*) ^< f]-^ >. 
The last texture contribution to these C' is, of course, independent of /, and 
proportional to < >. Hence we can expect a white noise type of spectrum 
(C' oc const) in this / region. In practice, the transition from white noise to scale 
invariance happens very quickly. We never observe a white noise regime, but only 
a suppression for the low /. 

This heuristic argument can be converted into an approximate formula by 
taking the toy-hat profile and carrying out explicitly the C' integration, starting 
from y =< yi >. Let us define li as 

16 



(43) 



Figure 8: Approximate spectra /(/ + 1)C' for various values of ps. 



(the /-scale of the last texture), and lis as 




47r 



(44) 



//.(//. + 1) 



(the /-scale of textures living at the last-scattering surface). Then 



C oc 




^1 — 



for / < h 
0'^ for / > h 



(45) 



giving mathematical expression to what we have said above. This leads to the 
approximate spectral formula 



with /i 0.5/ps and lis ~ for standard recombination, with ps = 2p^. 

The approximate spectra were plotted in Fig. 8 with an arbitrary vertical scale. 
Confronting this with Fig. 4, we see that (46) does provide a good qualitative 
description of the spectra. 

4 The two-point correlation function 

The two-point correlation function C{6) is a tool sadly closer to experiment than 
the C' spectrum (however see [14] for a more hopeful point of view). We will 
now show that C{6) is a very bad discriminator between inflationary and texture 
scenarios. C{6) is defined as 




for / < h 
for / > h 



(46) 




(47) 
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with 6 = (0,0'), (0 and 0' any two fixed directions in the sky). Making use of 
(17) this can be written as 

91 A- 1 

Ci9) = j:&^P,icos9). (48) 

The theoreticai C[9) can be estimated from the sky average of y"(^)^(^') with 
a fixed e = (OTO'): 

CoUO) = J ^^^(0)^(OXcos^-cos(0':0')) (49) 
which after some algebra leads to 

2/ -I- 1 

CobsiO) = Y.Ci——Piicos 9) (50) 
/ 47r 

with Ci defined as in (18). An expression for C{9) in texture scenarios can be 
found by replacing (22) in (48). It is curious that a simpler expression can be 
derived. Inverting the Legendre transform (15) one has 

91 A- 1 

W%z,y) = J2——W^\y)P,{z). (51) 
/ 47r 

If we define the Legendre-squared spot profile to be 

91 -A- 1 

W<'\z,y) = j:^W^''{y)P,{z) (52) 

then the correlation function is simply 

C{9) = —-^ ^dyN{y)W^('\co,9,y). (53) 
47r Jo 

The profile is the Legendre analogue of W^-kW^ (where -k is the convolution) 

in Fourier analysis. Depending on whether the profiles are defined in the 9 or 
in the Legendre spaces it may be easier to use (53) or to combine (48) and 
(22). If one wants to take into account the telescope beam filtering, one should 
multiply the C' in (48) by the square of the beam Legendre transform JF'^. In 
the numerics we will use an approximation to the COBE beam JF' = e~2('+2)^'^^ 
with (T = 27r(10V360°). 

We have also derived an expression for the cosmic variance in Cohs{9). In 
Gaussian theories ([4] with a factor of 2 correction) we have: 

'yl{Co,s{9)) = E^g(C0(^) Pi (cos 9) (54) 
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2- 



Figure 9: e(^) for Gaussian spots with (n = 1, = 0.25). 



with crQ[Ci) given by (36). In the texture scenarios not only does one have to 
replace (t'^{Ci) by (Tjxi^i) (^^ (39)), but also terms in P/P/' (with / ^ I') 
appear, reflecting the inter-/ correlations. A rather elaborate calculation, in the 
moulds of Sec. 3. 3, brings us to 



TX 



,,, 47r 47r 



Pi{cos 9) Pi, [cos 9) 



(55) 



with 



21 + 1 



+ 1 + 



J^^'dyN{y)W'^Hy)W''^'{y) 



< a' >V iltdy N{y)W'^'{y)){!tdy N{y)W''^'{y)) 



(56) 



where, again, we have neglected terms in l/Nf. The intensity of the non-Gaussian 
correction to the cosmic variance can be measured by means of the quantity 



TX 



;{Cobs{9)) 



(57) 



which depends on n, p^, and the profile. The dependence on n is trivial, so 
we define e{ps,9) = e{n,ps,9).n. The numerator in (57) can be rewritten as 
< (V|/«(2))2 showing that the variance in Cobs{9) is always larger in texture 
scenarios than in similar Gaussian theories. As before, we adopt the attitude 
that while this does reduce the predictability of the theory on C(^), it also signals 
non-Gaussian behaviour and the need of a non-Gaussian data-analysis approach 
to fully make out the predictions of the theory. In particular, on angular scales 
where e ^ 1, we know that it is worth studying the collapsed (to two points) 
n-point correlation function. In Fig. 9 we divided 9 G (0, tt) into 100 points 
and plotted e(^) for Gaussian spots with (n = 1, = 0.25). The fact that the 
low / multipoles contribute to C{9) for all 9 distributes non-Gaussianity nearly 
uniformly over 9. In Fig. 10 we analyse the various contributions to cr|.jf, sep- 



19 




Figure 10: cr^x (top line), and its Gaussian (bottom line), diagonal (top points), 
and ofF-diagonal (bottom points) in V;;/, contributions. 




Figure 11: C{6) for Gaussian spots with (n = 1, = 0.25) with its full cosmic 
variance error bars, confronted with its fitting n = 1.001 theory (top line at 
C(vr)). 

arating its Gaussian, diagonal and off-diagonal components in Vul We note that 
the inter-/ correlations responsible for the off-diagonal elements of Vu' can act 
so as to reduce the cosmic variance, an important fact which we shall make use 
of in [10]. Naturally, the subtle differences in the C'-spectra pointed out at the 
end of Sec. 3. 4 are completely drowned in C{6). As an example we have plotted 
in Fig. 11 C{6) for the above texture theory, and for its fitting tilted spectrum 
infiationary theory. As C{6) spreads the low / over all ^, the two theories come 
out completely confused. Note however that the normalization in the C' and in 
C{6) can be substantially different (in Fig. 11 we used a different normalization 
procedure than in Fig. 7). 
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5 The last texture and non-Gaussianity at low 



In the derivation of (22) we made use of statistics for unordered textures. However, 
a similar expression could have been obtained using ordered textures. In this case 
one should be wary for the W^' are not independent variables. Nevertheless for- 
mula (21) remains valid if one uses the marginal distribution functions of j/„ in 
computing < W^'^ > for each term in (21). The marginal distribution functions 
PniVn) are defined as 



PniVn) = J dyi...dyn-idyn+i...dyNtPi...Nt{yi,---,yNt) (58) 
and using (10) they can be found to be 



Pn{yn) = N{y^)e-^(y-^^^^^ 



(59) 



where M[y) = N[x)dx. Then: 



C = /"'^y.iV(y.)I^'-(y.)e-^(-)^p^ (60) 

47r ^^^Jo {n - 1)1 

and since Nt is very large we recover (22). The advantage of (60) is that it allows 
us to write 

C = (61) 

n 

with 

CL = ^ rdy. iV(y.)l^'-(y.)e--(-)^p^ (62) 
47r Jo [n — Ij! 

reporting the contribution of the n^^ texture to the multipole C'. In Figure 12 we 
plotted the C'* spectrum for four plausible scenarios: for a Gaussian profile with 
Ps = 0.05 and = 0.1 (n = 1 in both cases), and for a SSSS profile with n = 0.25 
and n = 0.1 (ps = 0.2 for both). Underneath the spectra we plotted the result of 
stopping the sum (61) at Nt from 1 to 10, thus obtaining the contribution of the 
first 10 textures to each C'. In all the cases considered we notice that the low 
/ multipoles are dominated by the contribution from the last texture. The / at 
which this stops being the case depends on the spot details (profile, n, p^), but 
/ = 2,3 always seem to be subject to this effect. Also the intensity of the last 
texture dominance at low / depends on the spot details, the effect being notably 
pronounced for Gaussian profiles with low ps. As we go up in /, the last texture 
becomes less prominent, but no other single texture replaces it as a dominant 
feature. We may ask how many textures are responsible for, say, 95% of a given 
C' for high /. Using the toy hat model and performing an approximate calculation 
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Gaussian profile with (n=l, p=0.05) 



Gaussian profile with (n=l, p=0.1) 
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SSSS profile with (n=0.25, p=0.2) 
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SSSS profile with (n=0.1, p=0.2) 
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Figure 12: Contribution of the last ten textures for a Gaussian profile, with 
Ps = 0.05 and = 0.1 (n = 1 in both cases), and for a SSSS profile with 
n = 0.25 and n = 0.1 (and ps = 0.2). 



similar to the one in Sec. 3. 5, we can show that this number is proportional to P. 
This implies that for high / the number of defects responsible for a given is 
proportional to /. 

These results have far reaching consequences. If one believes in fiippant ap- 
plications of the central limit theorem, then one can expect the for high / < 30 
to be independent Gaussian distributed random variables, as they add up a large 
number of independent contributions (proportional to /). We have now shown 
that this argument can certainly not be repeated for the low /. In fact, our remark 
in Sec. 3. 3 that cr'^xi^i) ^ ^Gaussi^i)} numerical results plotted in Fig. 

3, prove that the are distinctly non-Gaussian at low /. This also implies that 
these are dependent random variables [11]. We can now understand better 
why this is so. The non-Gaussian features in the for low / result from the 
fact that they are mostly due to one single lumpy object blatantly different from 
Gaussian noise. We may even expect the for low / to reproduce the morpho- 
logy, size, and other features of the last defect. Devising a data- analysis method 
capable of uncovering these features is the purpose of [10]. We should stress that 
both the central limit theorem and Fig. 3 strongly suggest, but still do not prove, 
that the are Gaussian for high / < 30. 
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Figure 13: The c? and c? function for a Gaussian (top), and a SSSS profile 



The non-Gaussian toy model studied in [11] showed in which way the low / 
multipoles would display non- Gaussian behaviour if they were fully due to the last 
texture. The contribution from n > 1 textures will soften this non-Gaussianity. 
Hence a plausible measure for the expected non-Gaussianity present in a low / 
multipole is a' = C\jC\ a quantity which depends on the profile, n, and p^. 
We have computed numerically c? and c? for a Gaussian and a SSSS profile. 
The results are plotted in Fig. 13. A better connection between a' and detailed 
features of non-Gaussianity measure distribution functions will be given in [10]. 



The focus in this paper was to derive for texture scenarios whatever one normally 
derives for Gaussian theories (C', C{9), and their cosmic variances). We met an- 
omalous behaviour (eg. formulae (39) and (55)) which (theoretically) indicates 
non-Gaussianity, but we have left to the sequels [10] and [15] the task of devel- 
oping a non-Gaussian data analysis method. On the whole, within the standard 
lore, the results obtained are surprisingly similar to infiationary scenarios. Still, 
we found a few novelties. We discovered that texture spectra are tilted, but 




6 What's new? 
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not uniformly. Fitting a tilted spectrum for / G (25, 30) (see Fig. 6) leaves a 
significant suppression of power at low / G (2,5) (see Fig. 7). We have also com- 
puted directly the cosmic variance error bars in C; and Co6s(^), and found them 
to be larger than in comparable Gaussian theories (see formulae (39) and (55) 
and Figs. 3 and fO). In fact, if the non-Gaussian correction is very large (which 
does happen for most of the parameter values suggested in the literature), cosmic 
variance drowns the otherwise I-cr suppression of power at low / mentioned above 
(see Fig. 7). However, one should not give way to despair. Large non-Gaussian 
corrections to cosmic variance only hint that we have not applied to the theory 
the right data analysis procedure. To give a flavour of [15] let us point out that 
as the Ci are dependent random variables in non-Gaussian theories ([If]), it may 
happen that the Ci spectrum shape seen by any observer is never the average C' 
spectrum shape. A more intelligent method to make predictions on C; spectra 
and to estimate global parameters in the C' spectrum is in order. It appears that 
even a proper study of the C' (or Ci) spectrum for textures requires non-Gaussian 
data analysis. 

The central result in this paper is undoubtedly the proo/ that texture low / 
multipoles are strongly non-Gaussian, and that the last texture is to be blamed 
for this (see Sec. 5). We showed how the low C' are mostly due to the last texture 
(slightly perturbed by the one after - see Fig. 12), and so we may expect the 
last texture's non-Gaussian features to be imparted on the low / multipoles. In 
[II] we proposed the use of m-structure and inter-/ measures to complement 
the CK m-structure measures act as multipole shape factors. Inter-/ measures 
correlate preferred directions in two multipoles. These two types of quantities 
are uniformly distributed in Gaussian theories. In topological defect scenarios we 
can expect low / shape factors to reflect the morphology of the last defect. Hence 
their distributions should be non uniform, peaking at different values for, say, 
texture and cosmic strings. Inter-/ measures should also display the correlations 
between the various low / which are all due to the same last defect. In remains 
to be seen ([10]) how much cosmic confusion there is between these signals and 
their Gaussian counterparts. 

On a lower key this paper was simply intended as an analytical model for 
texture CMBR skies with the greatest possible generality. We tried to leave as 
free parameters whatever simulations have failed to determine. The point is to 
find out what exactly must be decided in order to answer a particular question. 
We have concluded, for instance, that the C' spectrum shape is insensitive to n or 
a, and depends very little on the spots' profile, once the identifications (33) and 
(35) have been made. The controlling parameter is p^ which fixes the spectrum 
tilt at intermediate / (see Fig. 6). The spectrum normalization, on the other 
hand, depends on all the parameters of the model (see Sec. 3.1), which must be 
evaluated before any prediction for the texture symmetry breaking energy or 
bias parameter can be made. Interestingly enough, we found that the low / non- 
Gaussianity effects depend only on n, p^, and profile (Sec. 5). These are not free 
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parameters of the theory, so one ought to decide on this important issue purely 
by means of hard work. 

The paradigm proposed here cannot easily be adapted to cosmic string scen- 
arios. Neglecting non-circular spot modes is, in texture scenarios, only slightly 
inaccurate. The string counterpart to circular spots, line discontinuities cut-off at 
horizon distance (as in [f6]), is definitely a much cruder approximation. More im- 
portantly, neglecting texture spot-spot correlations introduces only a small error 
(in fact it would not be difficult to incorporate small correlations in our model). 
On the contrary, the various segments of the string network form a Brownian 
random walk. Segments linked to each other are therefore very highly correlated, 
making an analytical treatment cumbersome. Uncorrelated string segments (as 
in [f6]) are not very realistic. Overall, although we can boast a more realistic set 
of assumptions than [f6], this is mostly due to the fact that textures are far sim- 
pler than strings. This also justifies why we can go further than [f6]. We treated 
with greater rigour (no "ergodic" assumptions) a much broader set of models. 
We also managed to treat exactly (within the model) the problem of cosmic vari- 
ance. Finally, greater sophistication in the statistical machinery opened doors to 
the issue of non-Gaussian data analysis in defect scenarios (developed further in 
[10, 15]). 
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